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Quark propagator in Landau and Laplacian gauges with overlap fermions 
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The properties of the momentum space quark propagator in Landau gauge and Gribov copy free 
Laplacian gauge are studied for the overlap quark action in quenched lattice QCD. Numerical cal- 
culations are done on two lattices with different lattice spacing a and the same physical volume. 
We have calculated the nonperturbative wave function renormalization function Z{q) and the non- 
perturbative mass function M(p) for a variety of bare quark masses and perform a simple linear 
extrapolation to the chiral limit. We focus on the comparison of the behavior of Z{q) and M(p) in 
the chiral limit in the two gauge fixing schemes as well as the behavior on two lattices with different 
lattice spacing a. We find that the mass functions M(p) are very similar for the two gauges while the 
wave-function renormalization function Z{q) is more strongly infrared suppressed in the Laplacian 
gauge than in the Landau gauge on the finer lattice. For Laplacian gauge, it seems that the finite 
a error is large on the coarse lattice which has a lattice spacing a of about 0.124 fm. 

PACS numbers: 12.38.Gc, 11.15.Ha, 12.38.Aw, 14.65.-q 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) is widely accepted 
as the correct theory of the strong interaction and the 
quark propagator is one of its fundamental quantities. By 
studying the momentum-dependent quark mass function 
in the infrared region we can gain valuable insight into the 
mechanism of dynamical chiral symmetry breaking and 
the associated dynamical generation of mass. At high 
momenta, one can use the quark propagator to extract 
the running quark mass Ij. 

Lattice QCD provides a way to study the quark prop- 
agator nonperturbatively. There have been several lat- 
tice studies of the momentum space quark propaga- 
tor using different fermion 
actions. The usual gauge for these studies has been 
Landau gauge, because it is a (lattice) Lorentz covari- 
ant gauge that is easy to implement on the lattice, and 
the results from the lattice Landau gauge can be easily 
compared to studies that use different methods. Finite 
volume effects and discretization errors have been exten- 
sively explored in lattice Landau gauge 0,^3. Unfortu- 
nately, lattice Landau gauge suffers from the well-known 
problem of Gribov copies. Although the ambiguity orig- 
inally noticed by Gribov ^3] is not present on the lat- 
tice, since in practice one never samples from the same 
gauge orbit twice, the maximization procedure used for 
gauge fixing does not uniquely fix the gauge. In gen- 
eral, there are many local maxima for the algorithm to 
choose from, each one corresponding to a Gribov copy, 
and no local algorithm can choose the global maximum 
from among them. While various remedies have been 
proposed ,13, ,14], they are either unsatisfactory or com- 
putationally very intensive. For a recent discussion of the 



Gribov problem in lattice gauge theory, see Ref. |l5| . 

An alternative app roach is to work with the so-called 
Laplacian gauge [ig- This gauge is "Landau like" in 
the sense that it has similar smoothness and Lorentz 
invariance properties but it involves a non-local 

gauge fixing procedure that avoids lattice Gribov copies. 
The gluon propagator has already been studied in Lapla- 
cian gauge in Refs. 0, 0| and the improved staggered 
quark propa gato r in Laplacian gauge in Ref. 0|. It has 
been shown that Landau and Laplacian gauges be- 
come equivalent in the perturbative (high-momentum) 
regime and this has been confirmed by numerical stud- 
ies HIIllil. 

In this paper we study the overlap quark propagator 
in the Laplacian gauge and compare the results with the 
Landau gauge to explore the effects of selecting a gauge 
condition free of Gribov copies. Unlike Asqtad fermions, 
the overlap formalism provides a fermion action which 
is free of doublers and preserves an exact form of chiral 
symmetry on the lattice. The latter feature makes over- 
lap fermions the action of choice for studying dynamical 
chiral symmetry breaking near the chiral limit. 

We also compare Laplacian gauge results on two lat- 
tices with the same physical volume and different lattice 
spacings a to explore the finite a error. In this work, 
the C(a^) mean- field improved gauge action is used to 
generate the quenched gauge configurations. 



II. GAUGE FIXING 

We consider the quark propagator in Landau and 
Laplacian gauges. Landau gauge fixing is performed by 
enforcing the Lorentz gauge condition, X]u^m^m(^) ~ ^ 
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on a configuration by configuration basis. For the tadpole 
improved plaquette plus rectangle (Liischer-Weisz [2ll |') 
gauge action which we use in the current work, we use 
the 0(0^) improved gauge fixing scheme, this is achieved 
by maximizing the functional j23 | , 
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where Ti and J-2 are 

•^1 = jETr{^A^(^) + t^/'(^)} 



(1) 



and 



respectively, and uq is the usual plaquette measure of the 
mean link. In this Fourier accelerated, steepest- 

descents algorithm is used to find a local maximum. 
There are, in general, many local maxima and these are 
called lattice Gribov copies. This ambiguity in principle 
will remain a source of uncontrolled systematic error. 

Laplacian gauge fixing is a nonlinear gauge fixing 
that respects rotational invariance, has been seen to be 
smooth, yet is free of Gribov ambiguity. It is also com- 
putationally cheaper then Landau gauge fixing. There is, 
however, more than one way of obtaining such a gauge 
fixing in SU(N) lattice gauge theory. There are three im- 
plementations of Laplacian gauge fixing employed in the 
literature: 

1. 9^(1) gauge (QR decomposition), used by Alexan- 
drou et al. 

2. 9^(11) gauge, where the Laplacian gauge transfor- 
mation is projected onto SU(3) by maximising its 
trace [T^ . 

3. 9^(111) gauge (Polar decomposition), the original 
prescription described in Ref. [ifl] and tested in 
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All three versions reduce to the same gauge in SU(2). 
For a more detailed discussion, see Ref. |19j. For SU(3) 
staggered quarks, the study in Ref. ^ indicate that 9^(1) 
and 9'^(II) gauge give very similar results, and 9^(111) 
gauge is very noisy. In this work we will only use the 
9^(11) gauge. 



III. QUARK PROPAGATOR ON THE LATTICE 

In a covariant gauge in the continuum, the renormal- 
ized Euclidean space quark propagator has the form 
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(2) 



where C is the renormalization point. The renormaliza- 
tion point boundary conditions arc chosen to be 



(3) 



where m(C^) is the renornialized quark mass at the renor- 
malization point. The functions A{('^;p'^) and B{('^;p'^), 
or alternatively Z{('^;p'^) and M{p'^), contain all of the 
nonperturbative information of the quark propagator. 
Note that M^p^) is renormalization point independent, 
all of the renormalization-point dependence is carried by 
Z{e;p^). 

When all interactions for the quarks are turned off, 
i.e., when the giuon field vanishes (or the links are set to 
one), the quark propagator has its tree- level form 



S^'\p) = 
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where rrfl is the bare quark mass. When the interactions 
with the gluon field are turned on we have 

S(^){p) ^ 5b"'=(a;p) = Z2{e:a)S{e;p), (5) 

where a is the regularization parameter - in this case, 
the lattice spacing - and Z2(C^;a) is the quark wave- 
function renormalization constant chosen so as to ensure 
■^(C^;P^)lp2=c^ ~ ^- simplicity of notation we sup- 
press the a-dependence of the bare quantities. 

On the lattice we expect the bare quark propagators, 
in momentum space, to have a similar form as in the 
continuum, except that the 0(4) invariance is replaced 
by a 4-dimensional hypercubic symmetry on an isotropic 
lattice. Hence, the inverse lattice bare quark propagator 
takes the general form 
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Y.^Mi)j +B{p). 
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With the periodic boundary conditions in the spatial di- 
rections and anti-periodic in the time direction, the dis- 
crete lattice momenta will be 

where rii — 1, .., Ni and nt — 1, .., Nt, Ni and Nt are the 
lattice extent in spatial and temporal direction respec- 
tively. 

The overlap fermion formalism |25| realizes an ex- 
act chiral symmetry on the lattice and is automatically 
0{a) improved. The massive overlap operator can be 
written as l2a 
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where Hw{x,y) = 75-Di«(a;, y) is the Hermitian Wilson- 
Dirac operator, e{Hw) = H^j ^JlT^ is the matrix sign 
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function, and the dimensionless quark mass parameter rj 
is 



,0 



2m„ 



(9) 



where m'^ is the bare quark mass and is the Wilson 
quark mass which, in the free case, must be in the range 
< rriuj < 2. The bare quark propagator in coordinate 
space is given by the equation 



(10) 



where 



2m, 



-D~^{rf) and 



1 — rj 



(11) 



When all the interactions are turned off, the inverse 
bare lattice quark propagator becomes the tree-level ver- 
sion of Eq. lO 

{S^'^)-\p)^^{Y,C^^\p)^^ +B^'Hp). (12) 



We calculate S'*-"-' (p) directly by setting the links to unity 
in the coordinate space, doing the matrix inversion and 
then taking its Fourier transform. It is then possible to 
identify the appropriate kinematic lattice momentum q 
directly from the definition 



(13) 



The form of qfj,{pfj,) is shown and its analytic form given 
in Ref. 9] . Having identified the appropriate kinematical 
lattice momentum q, we can now define the bare lattice 
propagator as 



S'""'''ip) = 



H + M{p) ■ 



(14) 



This ensures that the free lattice propagator is identical 
to the free continuum propagator. Due to asymptotic 
freedom the lattice propagator will also take the contin- 
uum form at large momenta. In the gauge sector, this 
type of a naly sis dramatically improves the gluon propa- 
gator miiilli. 

The two Lorentz invariants can then obtained by 



M{j>)^^Tt{S-\p)}. 



(15) 
(16) 



This means that Z{p) is directly dependent on our choice 
of momentum, g, while M{p) is not. 



IV. NUMERICAL RESULTS 
A. Simulation parameters 

In this paper we work on two lattices with different 
lattice spacing, a, and very similar physical volumes. 
The gauge configurations are created using a tadpole 
improved plaquette plus rectangle (Liischer-Weisz pl|) 
gauge action through the pseudo-heat-bath algorithm. 
For each lattice size, 50 configurations are used. Lat- 
tice parameters are summarized in Table ^ The lattice 
spacing a is determined from the static quark potential 
with a string tension ^ = 440 MeV [soj . 

Landau gauge fixing to the gauge configuration 
was done using a Conjugate Gradient Fourier Ac- 
celeration ^3 algorithm with an accuracy of 6* = 



Y.\dM^)\ < 10- 



The improved gauge-fixing 



scheme was used to minimize gauge-fixing discretization 
errors ^251. For the Laplacian gauge fixing, we only use 
the (II) gauge J_9] . In this case we construct the gauge 
transformation by projecting M{x) constructed from the 
three lowest lying eigenmodes, onto SU(3) by means of 
trace maximisation. Effectively, we maximise the trace 
of G{x)M{xy by iteration over Cabibbo-Marinari SU(2) 
subgroups. 

Our numerical calculation begins with an evaluation of 
the inverse oi D{rf} with the unfixed gauge configurations, 
where D{rf) is defined in Eq. 0. We approximate the 
matrix sign function e{Hm) by the 14th order Zolotarev 
approximation j32|J- We then calculate Eq. H10() for each 
configuration and rotate it to Landau or Laplacian gauge 
by using the corresponding gauge transformation matri- 
ces {Gi{x)}. Afterward we take the ensemble average to 
obtain 3^'^'^° {x , y) . The discrete Fourier transformation 
is then applied to S^'^'''^ (x , y) and the momentum-space 
bare quark propagator, S^'^'^°{p) is obtained finally. 

We use the mean-field improved Wilson action in the 
overlap fermion kernel. The value k = 0.19163 is used 
in the Wilson action, which provides m^ja = 1.391 for 
the regulator mass in the interacting case Q . We calcu- 
late the overlap quark propagator for ten quark masses 
on each ensemble by using a shifted Conjugate Gradient 
solver. For the two lattices considered here, the quark 
mass parameter 77 was adjusted to make the tree level 
bare quark mass in physical units, the same on the two 
lattices. For example, we choose jj, — 0.018, 0.021, 0.024, 
0.030, 0.036, 0.045, 0.060, 0.075, 0.090, and 0.105 on 
ensemble 1, i.e., the 16'^ x 32 lattice with a = 0.093 
fm. This corresponds to bare masses in physical units 
of m° = 2^TO„ = 106, 124, 142, 177, 212, 266, 354, 442, 
531, and 620 MeV respectively. 

The results of lattice 2 (12'^ x 24) in Landau gauge were 
presented in detail in Ref. , and the results of lattice 1 
(16'^x32) in Landau gauge were also reported in Ref. [ll| . 
Here we will focus on the comparison of the results of two 
lattice gauge fixing schemes, i.e., the Landau gauge and 
the Gribov copy free Laplacian gauge, to probe the be- 
havior of the overlap fermion propagator with different 
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TABLE I: Lattice parameters. 



Action 


Volume 


A'rhcrm 


A^'samp 


P 


a (fm) 




Physical Volume (fm*) 


Improved 


16^ X 32 


5000 


500 


4.80 


0.093 


0.89650 


1.5^ X 3.00 


Improved 


12^ X 24 


5000 


500 


4.60 


0.124 


0.88888 


1.5^ X 3.00 



GeV 



5 6 7 
q GeV 



10 11 12 



FIG. 1: (Color online). The functions M{p) and Z^^\q) = 
■^(C^i q) renormalized at ^ = 5.31 GeV (in the q scale ) for all 
ten quark masses in Laplacian gauge on the 16'^ x 32 lattice. 
The mass function M{p) is plotted versus the discrete mo- 
mentum p defined in Eq. Q, p = v'X] P^u c>"^6r the interval 
[0,5] GeV, and Z^^\q) is plotted against the kinematic mo- 
mentum q defined in Eq. II13II . q = v'X] 9?' over the interval 
[0,12] GeV. The data correspond to bare quark masses (from 
bottom to top) ^Ji = 0.018, 0.021, 0.024, 0.030, 0.036, 0.045, 
0.060, 0.075, 0.090, and 0.105, which in physical units corre- 
spond to mP = 2riTn^ ~ 106, 124, 142, 177, 212, 266, 354, 
442, 531, and 620 MeV respectively 



B. Laplacian gauge 

In Fig. n we show the results for all ten masses for 
both the mass and wave-function renormalization func- 
tions, M{p) and = Z{(^;q) respectively. As was 
shown in Ref. the continuum limit is more rapidly 
approached when the mass function is plotted against 
the discrete lattice momentum p, while the wave-function 
renormalization function Z^^^(q) is plotted against the 
kinematic momentum q. The renormalization point in 
Fig. □ for ^("'(q) has been chosen to be ( = 5.31 GeV 
in the q-scale. 

In the plots of M{p), the data is ordered as one would 
expect by the values for bare quark mass mP, i.e., the 
larger the bare quark mass mP, the higher the M{p) 
curve. For large momenta, the function Z^^\q) demon- 
strates little mass dependence. Deviation of Z^^^ (q) from 
its asymptotic value of 1 is a sign of dynamical symmetry 
breaking, so we expect the infrared suppression to van- 
ish in the limit of an infinitely heavy quark. In the figure 
for Z^^\q), the smaller the bare mass, the more pro- 
nounced is the dip at low momenta. Similarly, at small 
bare masses M{q) falls off more rapidly with increasing 
momenta, which is understood from the fact that a larger 
proportion of the infrared mass is due to dynamical chiral 
symmetry breaking at small bare quark masses. These 
results are much the same as in Ref. which display 
the data on the same lattices in Landau gauge. This 
qualitative behavior is also consistent with what is seen 
in Dyson-Schwinger based QCD models [sslls^. 

In Fig. 13 we plot the data after a linear chiral extrap- 
olation for both functions M{p) and Z^^\q) = Z{(;q) 
in Laplacian gauge. The mass function M{p) is shown 
against p and while the wave function renormalization 
function Z^^\q) is shown against q with the renormal- 
ization point chosen as at 5.31 GeV in the q scale. We 
see that both M{p) and Z^^^ [q) deviate strongly from the 
tree-level behavior, which are M{p) = mP and Z^^\q) = 
1. In particular, as in earlier studies of the Landau gauge 
quark propagatorH H IE H 113, 

we find a clear signal 
of dynamical mass generation and a significant infrared 
suppression of the Z{(; q) function. 



gauge fixings, and the effect of Gribov copies. Before we 
make the comparison, we first briefly present some data 
in Laplacian gauge on the 16'^ x 32 lattice, our fine lattice. 
All data has been cylinder cut . Statistical uncertain- 
ties are estimated via a second-order, single-elimination 
jackknife. 



C. Gauge fixing comparison 

Next we present the results in the two gauge fixing 
schemes for comparison. All data have been cylinder 
cut [23|- First we give the results on the 16^ x 32 lat- 
tice. 
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FIG. 2: The plot of the functions M{p) and Z'^^^q) = Z{C\ q) 
after a hnear extrapolation to the chiral limit. The mass func- 
tions M{jp) is plotted against the discrete momentum p in the 
upper part of the figure and Z''^\q) with the renormaliza- 
tion point C, =5.31 GeV (in the g-scale) is plotted against the 
kinematic momentum q in the lower part of the figure. 



Fig-Elreports results for the mass and renormalization 
functions at our lightest bare quark mass of wP = 106 
MeV. These results may be compared with the Asqtad 
results of Ref. where Fig. 9 compares results of var- 
ious gauge fixing schemes for the renormalization func- 
tion. There, Landau gauge results are seen to lie sig- 
nificantly higher than the Laplacian gauge results in the 
infrared. However, with overlap fermions, Fig.|21indicates 
the Landau gauge results lie much closer to the Lapla- 
cian gauge results. Given the improved chiral properties 
of the overlap operator, the present results should pro- 
vide a better indication of the continuum limit behavior. 
In either case, the same qualitative behavior of Landau 
Gauge sitting above Laplacian gauge in the infrared is 
observed. 

The mass function of Fig.|3|reveals an approximate in- 
variance on the selection of Landau or Laplacian gauges. 
Figs. 12 and 13 of Ref. 6] indicate that the mass func- 
tion of Asqtad fermions is also insensitive to the choice of 
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FIG. 3: (Color online). The comparison of two gauge fixing 
results at finite bare quark mass (m° = 106 MeV) on the 
fine 16^ X 32 lattice with a = 0.093 fm. Z'-^^q) is renor- 
malized to one at the renormalization point ( =5.31 GeV 
(in the q-scale). For the mass function M(p), Landau gauge 
and Laplacian gauge are very similar, while the wave func- 
tion renormalization functions Z^-^^ (q) differ in the infrared 
region. 



Landau gauge or the Gribov-copy free Laplacian gauge. 

We now proceed to compare the data in the chiral limit. 
Fig.^lshows the comparison of the mass function, M{p), 
and the wave function renormalization function, Z^-^'^{q), 
in Landau gauge and Laplacian gauge. We see that they 
give similar performance in terms of rotational symme- 
try and statistical noise. Looking more closely, we can see 
that Landau gauge gives a slightly cleaner signal at this 
lattice spacing. We also note that at very large momenta, 
the two gauge fixing schemes give similar results as ex- 
pected. Although Laplacian gauge is a non-local gauge 
fixing scheme and difficult to understand perturbatively, 
it is equivalent to Landau gauge in the asymptotic re- 
gion [23 • In the infrared region, the mass function, M{p), 
in the two gauges are very similar. For the mass func- 
tions there is a hint that the data in Laplacian gauge are 
a little higher than for the Landau gauge, although they 
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FIG. 4: (Color online). The comparison of two gauge fixing 
results after a linear extrapolation to the chiral limit on the 
fine lattice, i.e., 16^ x 32 at a = 0.093 fm. Z'-^^{q) is renor- 
malized to one at the renormalization point =5.31 GeV 
(in the g-scale). For the mass function M{p), Landau gauge 
and Laplacian gauge are very similar, while the wave function 
renormalization functions Z'^' (g) are similar in the large mo- 
mentum region but differ in the infrared region. 
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FIG. 5: (Color online). The comparison of two gauge fixing 
results in the chiral limit on the coarse lattice, i.e., 12'^ x 24 
with a = 0.124 fm. The small gauge dependence of the in- 
frared behavior of the Z-function is similar to that on the fine 
lattice in Fig. 2] The infrared mass functions, M(p), appear 
different on this coarse lattice whereas they were similar on 
the fine lattice. This suggests larger C'(a^) errors in Laplacian 
gauge. 



agree within statistical errors. With greater statistics we 
may resolve a small difference. For the renormalization 
function Z'^^^q), there are systematic differences in the 
infrared region. The Z^^\q) is more strongly infrared 
suppressed in the Laplacian gauge than in the Landau 
gauge. That is consistent with what was seen in the case 
of the Asqtad quark action 6] when comparing these two 
gauges. 

Now we give the results from the coarse lattice. Fig.|Sl 
shows the comparison of the mass function M{p) and 
the wave function renormalization function ^'•^•'(g) in 
the chiral limit in Landau gauge and Laplacian gauge on 
the 12^ X 24 lattice with a = 0.124 fm. As in the case of 
the 16'^ X 32 lattice, at very large momenta, the two gauge 
fixing schemes give similar results. For the renormaliza- 
tion function Z^^'^{q), the situation is very similar to the 
case of the 1& x 32 lattice, i.e., at very large momenta. 



the two gauge fixing schemes give similar results, while in 
the infrared region, Z^^^ (q) is more strongly suppressed 
in Laplacian gauge than in Landau gauge. For the mass 
function, M(p), the situation is different to that seen on 
the fine lattice. In the infrared region, the data in Lapla- 
cian gauge sit higher than in Landau gauge. While this 
is consistent with the infrared behavior of the renormal- 
ization function, Z^^\q), which is more suppressed in 
Laplacian gauge than in Landau gauge, there is no sim- 
ilar signal for our fine lattice 16"^ x 32. In that case (see 
Fig. O , the data in Laplacian gauge agree with that in 
Landau gauge within error bars, although there is a hint 
that the data in Laplacian gauge are a little higher than 
in Landau gauge in the infrared region. This infrared be- 
havior of the mass function in Laplacian gauge is likely 
caused by the finite lattice spacing errors, i.e., on fine 
enough lattices we expect Landau and Laplacian gauge 
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well at large momenta, but there is a substantial differ- 
ence in the infrared region. A similar comparison was 
made between these two lattices in Landau gauge [llj . 
In that case, both the renormalization function, Z^^''>{q), 
and the mass function, M{p), agree well on the two lat- 
tices. This indicates that in Landau gauge, the finite a 
errors are small even for our coarse lattice with lattice 
spacing a — 0.124 fm. But in Laplacian gauge, the finite 
a errors are not negligible for our coarse lattice 12^ x 24. 
A likely explanation for this is the fact that our imple- 
mentation of the Laplacian gauge is not C(a^) improved. 
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FIG. 6: (Color online). The comparison of Laplacian gauge 
results in the chiral limit on two lattices, i.e., 12'^ x 24 with a 
= 0.124 fm and 16^ x 32 with a = 0.093 fm. The renormahza- 
tion point for Z'^'(q) = Z{(; q) is chosen to be C =5.31 GeV 
(in the g-scale). For the mass functions M{p), the value on 
12^ X 24 lattice is higher than that on the 16^ x 32 lattice in 
the infrared region, while the wave function renormalization 
function Z^^^ (g) on two lattices agree within error bars. This 
suggests that the finite a errors are not small on the coarse 
lattice in Laplacian gauge. 



V. SUMMARY AND OUTLOOK 

The momentum-space quark propagator has been 
studied in Landau gauge as well as the Gribov copy free 
Laplacian gauge on two lattices with the same physical 
volume but with different lattice spacings a . We calcu- 
lated the nonperturbative momentum-dependent wave- 
bmction renormalization, Z{q), and the nonperturbative 
mass function, M{p), for a variety of bare quark masses. 
We also performed a simple linear extrapolation to the 
chiral limit. 

At very large momenta the two gauge-fixing schemes 
give similar results as expected. Laplacian gauge is equiv- 
alent to the Landau gauge in the asymptotic region. In 
the infrared region, the mass function, M{p), in the two 
gauges are very similar on the fine lattice, but differ on 
the coarse lattice. The present Laplacian gauge fixing is 
not C(a^) improved. For our fine lattice in the infrared 
region, the mass function, M{p), agrees within statistical 
errors in the two gauge fixings. However, there is a hint 
that the data in Laplacian gauge may be a little higher 
than in Landau gauge, which would be consistent with 
the behavior of the renormalization function, Z^^^{q). 
For the renormalization function, Z(^-'(q), there are sys- 
tematic differences in the infrared region. The renormal- 
ization function is more strongly infrared suppressed in 
the Laplacian gauge than in the Landau gauge. 



results for M{p) to be very similar. 

Finally, we present the data for the two lattices in 
Laplacian gauge to further explore possible finite a er- 
rors. Fig. shows the comparison of the mass function, 
M{p), and the wave function renormalization function 
2'(f^)(g) in the chiral limit in Laplacian gauge on the 
12'^ X 24 and the Idi^ x 32 lattices. For the renormal- 
ization function, Z^^^(g), results from the two lattices 
have small differences, but agree with each other within 
errors. For the mass function, M(p), the results agree 
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